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WEAK LOCALLY HOMOGENEOUS TURBULENCE IN 
IDEALIZED FLOW THROUGH A CONE 
by Robert G. Deissler 
Lewis Research Center 

SUMMARY 

Turbulence in idealized flow through a cone is analyzed by using generalized two- 
point correlation equations that are based on the Navier-Stokes and continuity equations. 
Viscosity and the effects of mean strain are retained in the analysis, but the turbulence 
is assumed to be weak enough to neglect the triple correlation terms in the equations in 
comparison with the other terms. Components of the turbulent velocity and vorticity 
variances, as well as the spectra of those quantities, are calculated. A mean-strain 
energy-transfer term in the spectral equation which transfers energy components 
between wave numbers is also considered. By using the momentum-heat-transfer anal- 
ogy, the results are related to heat transfer between the fluid and the cone wall. 

INTRODUCTION 

When fluid flows through a section of a cone, the fluid elements are distorted because 
of the changing cross sectional area of the flow. If turbulence is initially present, it is 
modified by this distortion of the fluid. The interaction between weak turbulence and 
an idealized distorting mean flow in a cone is studied in this report. Flow in a cone is 
of particular interest because wind tunnels and rocket nozzles frequently contain conical 
sections. 

The effect of an irrotational distortion (no shear) on turbulence has been studied in 
references 1 to 4. In these studies the effects of viscosity were neglected and the dis- 
tortion was assumed to occur so rapidly that the turbulent velocities have a negligible 
effect on the motion during distortion. The present work is more closely related to that 
of Pearson (ref. 5). In his work the effects of viscosity are included*and the requirement 
that the distortion be rapid is not imposed. If the distortion is not rapid, however, the 
turbulence must be weak enough in comparison with other terms in the equations to 



neglect terms containing triple correlations. Pearson's analysis assumes that the normal 
velocity gradients are uniform and that the turbulence is homogeneous. 

The present analysis of turbulence in incompressible idealized flow through a cone 
uses the generalized two-point correlation equations that are based on the Navier-Stokes 
and continuity equations. The normal velocity gradients 9Uj/9xj, 9 U 2 / 9 X 2 and 
9Ug/9xg are allowed to vary axially but not transversely. (Symbols are defined in the 
appendix. ) The turbulence is assumed to be homogeneous in the transverse direction 
but only locally homogeneous in the axial direction. That is, the variation in the intensity 
of the turbulence over a correlation (or mixing) lei^h in the axial direction is assumed 
to be small. The mean axial velocity is assumed uniform over a cross section, and 
mean shear stresses are assumed absent. The turbulence is initially isotropic but is 
allowed to become anisotropic under the distorting influence of the mean flow. Compo- 
nents of the turbulent velocity and vorticity variances are calculated, as well as compo- 
nents of the spectra of those quantities. A mean-strain energy-transfer term in the 
spectral equation which transfers energy components between wave numbers is also 
considered. By using the momentum- heat-transfer analogy, the results are related to 
heat transfer between the fluid and the cone wall. The basic equations required in the 
analysis will be considered in the next section. 

BASIC EQUATIONS 

General two-point correlation equations for turbulence in an incompressible fluid 
with mean velocity gradients are obtained in reference 6 from the Navier-Stokes and 
continuity equations as follows: 
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The vector configuration for correlations between fluctuatit^ quantities at points P and 
P* is shown in figure 1. The quantities u^ and ut are fluctuating velocity components 
at P and P , and Uj are mean velocity components, is a space coordinate, t 
is the time, p is the density, v is the kinematic viscosity, and p is the instantaneous 
pressure. Bars over terms designate correlations or averaged quantities. The sub- 
scripts can take on the values 1, 2, or 3, and a repeated subscript in a term indicates a 
summation. 

For locally homogeneous turbulence the turbulence is considered homogeneous over 

a correlation length, or the scale of the inhomogeneity 
is much greater than the scale of the turbulence. Thus, 
a quantity such as 02u.u'j/ 9(xjj.)jj^0(xj^)jjj^ in equation (1) 

will be negligible compared with 0^'ij^/9rjj0rj^. (A 
calculation for axially decaying turbulence without mean 
velocity gradients (ref. 6, fig. 3) shows that this is a 
good approximation except very close to the virtual 
origin of the turbulence. ) In general, for locally homo- 
geneous turbulence. 



Ffgure 1. - Vector configuration for two-point 
correlation equations. 
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Also, for that type of turbulence the mean velocity will vary linearly over distances for 
which correlations are appreciable so that 
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Finally, in order to make the set of equations determinate, the turbulence is assumed 
weak enough to neglect terms containing triple correlations. It should be noted that the 
turbulence in a flow with large velocity gradients may not have to be as weak as that in 
a flow without velocity gradients. The terms containing those gradients may be large 
compared with triple correlation terms, even if the turbulence is moderately stroi^. 
Equations (1), (2), and (3) become, for the steady state at a fixed point, 
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Equations (4), (5), and (6) are the correlation equations for steady-state locally 
homogeneous turbulence with mean velocity gradients. The equations can be converted 
to spectral form by introducing the usual three-dimensional Fourier transforms defined 
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as follows: 
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where is the Kronecker delta. Equation (10) can be obtained by differentiating 
equation (7) with respect to rj^, writing the inverse transform, and then dlfferentiatii^ 
with respect to . Taking the Fourier transforms of equations (4) to (6) results in 


^kj \ *^ik 


9(UJ 


k'm 


8(%) 


k^m 


9(Xk)m 


‘i'm 


d(0 . . 

dK„ 


+ ®£k‘^ij +^^k^m 


3(Xk) 


(P 


ij 


m 


^ (i/CjXl - i/c.\.) - 2vk\^^ (11) 


- i iK.Xl = 2 -LJ cp 




( 12 ) 


Combining equations (11) and (12) and noting that = 3Uj^/3Xj^ = 0 by 

continuity result in 
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Equation (13) is the spectral equation for steady-state weak locally homogeneous 
turbulence with mean velocity gradients. 

Consider next the case where the mean strain is irrotational, that is, where the 
shearing components of the mean velocity gradient are zero. Let 
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Then equation (13) becomes, for irrotational strain. 
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where (f ) is not a tensor subscript. For axisymmetric strain at each point in a cross 
section, as occurs in uniform flow through a cone, s^ 2 ) = ®( 3 )> continuity of the 

mean flow. 
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The subscript 1 refers to the direction of an axis of symmetry for the turbulence. Since 
S(i) = f(xj^), integration of equation (15) gives U 2 = - (l/ 2 )s^j^^X 2 and U 3 = - (l/2)s^jjXg, 
or, for a circular cross section, the radial velocity Uj.= - (l/2)s^j^^r. In addition, it is 
assumed that the turbulence is homogeneous over a cross section of the flow and that the 
turbulence changes only in the axial or Xj direction, so that 
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To simplify the notation, let (Uj)^^ = U, (Xj)^^^ = x, and s^j^^ = s in the remainder of 
this report. Then for equation (14) becomes 
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where use was made of the continuity relation in the form ^ 2 * 5 ^ 12 '^3*^13 “ " '^l^ll* 
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Next determine U and s as functions of x for flow through a cone of arbitrary 
cross section. With the aid of the diagram in figure 2 these are obtained as 
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Figure 2. - Section of cone considered in analysis. Flow cross 
section is arbitrary; a can be positive or negative. 
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Then equations (17) and (18) become 


d(p 


11 


X — + 2k^ — ^ - K 


9^^!! 


ax 


dK-i 


BKc 


a«. 


}- 12 ^ + 4 

k2 


xVjcpii 

Vo 


( 22 ) 


and 


X 2«, ^ ^ - 2/l*£J2!!l“xVWu.6^„ 


ax 


ax. 


a«r 


dKr 


^0^0 


(23) 


Equations (22) and (23) apply to either diverging or converging flow through a cone, 
depending on whether a is positive or negative. 


SOLUTIONS OF SPECTRAL EQUATIONS 

Before equations (22) and (23) can be solved in an initial value problem, the turbu- 
lence must, of course, be specified at an initial position. It is assumed that the turbu- 
lence is isotropic at x = Xq, the virtual origin of the turbulence (or at X = Xq = r^/tan a, 
by eq. (21)), and that 


(<Pij)0 = - «i/Cj) (24) 

12 JT^ 

where Jq is a constant (ref. 6, eq. (43)). 

It appears that the use of equation (24) is reasonable because it is the simplest con- 
dition that gives results that, at all values of x, reduce to those for isotropic turbulence 
as the mean strain goes to zero. The use of that initial condition implies that Pearson' s 
parameter t = vk^s approaches <», where Kq is a characteristic initial wave number 
of the turbulence (ref. 5). Thus, the present results should be applicable for large 
kinematic viscosity, small initial turbulence scale, or small strain rate. The case 
T — oo was not considered by Pearson. 

Equations (22) and (23) are first-order partial differential equations in four inde- 
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pendent variables and can be solved by methods similar to those given in reference 7. 
Solutions of these equations in dimensionless form, subject to the initial condition given 
in equation (24), are 
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It is of interest that and are functions only of (x - Xq) tan oi/tq and k*. 

Because of axial symmetry it is not necessary to obtain equations for <^22 ‘P33‘ 

In order to integrate over wave-number space, it is convenient to introduce spher- 
ical coordinates as follows: 
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Equation (7) then becomes, for r = 0, 
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The quantity xf/^ is a function only of the magnitude of the wave number k and 
represents an energy spectrum tensor which has been integrated over all directions in 

wave-number space. 

The velocity variances u^, u^u., Ug, and Ug are calculated from equations (25) to 
(30) as 
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and X* is given by equation (28). 

An important physical quantity related to the turbulence is the vorticity tensor oTTa^. 
The vorticity spectrum tensor is given in reference 8 as 


^ii = (42) 

By analogy with equation (33) a directionally integrated vorticity spectrum tensor can be 
defined as 
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Then the vorticity tensor is given by 


(43) 


11 



(44) 


W.CO. 
1 ] 



d/C 


One other quantity of considerable interest is the transfer term in equation (11) 
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This term is, of course, nonzero only in the presence of mean strain. The term was 
previously discussed in references 6 and 9 for the case of a shear flow. That it can be 
interpreted as a transfer term can be seen from equation (10), where, if r= 0, 



(46) 


That is, when the term is integrated over all wave-number space it gives zero contribu- 
tion to the rate of change of u.u.. (The quantity (Ujj)^^ ^ 

interpreted as a rate of change. ) The term /3j^j can, however, transfer energy between 
wave numbers. Evidently the transfer takes place by the stretching or compressing of 
the vortex lines associated with the turbulence by the mean strain. This transfer is 
similar to that produced by triple correlations, with the exception that in that case the 
stretching or compressing of the vortex lines is accomplished by the action of the turbu- 
lence on itself, rather than by a mean strain. For the present case of axially symmetric 
irrotational strain, equation (45) becomes, for i = j = 1, 
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where is obtained from equation (25). As in the case of (eq. (33)), can be 
integrated over all directions in wave-number space to give 
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Computed velocity variances, vorticities, and spectra will be discussed in the next 
section. 
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Figure 3. - Variance of turbulent velocity components for flow through cone. 


RESULTS AND DISCUSSION 

Velocity and Vorticity Variances 

2 2 2 

Turbulent velocity variances Uj^, Ug, and Ug, calculated using equations (34), (35), 

is divided 
ed is the 

curve obtained by solvii^ equation (22) as though the effects of strain were absent by 
omitting the transfer, pressure, and production terms (terms 2 to 6). The resulting 
equation is 


and (36), are plotted in dimensionless form in figure 3. In the ordinate, u. 
by X’ in order to give quantities which are independent of x. Also incl 


Uq tan a 
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Negative values of x* , the abscissa, correspond to a converging flow and positive values 
to a divergii^ flow. The virtual origin of the turbulence (x’ = 0) is the point at which the 
energy of the turbulence would be infinite. 

In a convergir^ flow, velocity fluctuations decrease from the virtual origin toward 
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the apex of the cone because of viscous dissipation. The distorting influence of the cone 
causes the longitudinal component to decrease more rapidly and the transverse compo- 
nents to decrease less rapidly than they would without the effects of strain. As the apex 
of the cone is approached, the longitudinal component continues to decrease rapidly 
toward zero. The transverse components, on the other hand, begin to increase as the 
effects of strain become greater than the effects of viscous dissipation. From equa- 
tion (20) it can be seen that the strain rate increases rapidly as the apex is approached. 

At the apex, where the mean velocity and strain rate approach infinity, the transverse 
velocity fluctuations would also become infinite. In practice, of course, the tip of the 
cone would have to be removed in order to allow flow. It is of interest that this increase 
in transverse turbulent velocity component with mean velocity or contraction ratio at 
large velocity ratios was not observed in reference 5, where the velocity gradients were 
uniform. 

For a diverging flow in a cone (positive x') the effects of strain are opposite to those 
for a converging flow. The longitudinal component decreases less rapidly than it would 
for no strain. Although the effect of strain is to increase the longitudinal component, 
that component continues to decrease as x* increases because, as shown in equation (20), 
the strain rate decreases with x* . That is, the effect of viscous decay in this case is 
greater than the effect of strain. 

Trends similar to those shown in figure 3 for converging flow have been observed 
experimentally (refs. 10 and 11). A comparison between the present analysis and 
experimental results from reference 10 for a 4-to-l contraction is given in fiigure 4. The 
subscript a refers to conditions upstream of the contraction at the point where the mean 
velocity begins to vary. (It should be noted that the contraction in the experiment was 
evidently not a true cone. ) Values of U/U_ for the analytical curves were obtained by 

A 



Figure 4. - Present analysis compared with experiment for converging flow. 
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Figure 5. - Variance of turbulent vorticity components for flow through 
cone. 


assuming that the value of U/U„ at the minimum point in the analytical curve for the 

a 

transverse component (i = 2) corresponds to the value of U/U^ at the minimum point in 
the experimental curve for the same component. The overall change in the turbulent 
components produced by the contraction appears to be given reasonably well by the 
analysis, but the minimum in the experimental curve for i = 2 is sharper than that in 
the corresponding analytical curve. 

To explain the trends shown in figure 3, turbulent vorticity variances are plotted in 
figure 5. The dashed curves in the figure, for no effects of strain, were obtained from 
the equation 
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In general, the trends for the vorticity components are opposite to those for the veloc- 
ities. For a converging flow, the longitudinal vortices are stretched and thus strengh- 
ened by the accelerating flow and the transverse vortices are shortened along their axes 
and are thus weakened. The straining action also tends to turn the axes of the oblique 
vortices and to aline them in the direction of the cone axis. Both of these effects tend 
to increase the loi^itudinal vorticity component and decrease the transverse component. 
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Figure 6. - Plot of u\jul and against velocity ratio. 

as shown in figure 5. The opposite trends occur for a divergii^ flow; the strain in that 
case tends to decrease the longitudinal vorticity component and to increase the transverse 
component. 

The relation of the vorticity to the fluctuating velocity components for a contracting 
flow has been pointed out by Prandtl (ref. 1) and Taylor (ref. 2), When the vortex axes 
lie predominately in the direction of the mean flow, the longitudinal velocity fluctuations 
are small, whereas the transverse velocity fluctuations can be large. The opposite 
effects occur for a diverging flow. There is, however, an important difference between 

2/2 2/2 

the two cases. This is illustrated in figure 6, where ^ 1/^2 plotted 

against U/Uq. (The relation between this abscissa and the one in the preceding plots is 

2/2 2/2 

given by eq. (19). ) Although uWUg goes to zero as w ^ /wg goes to infinity for a con- 

2/2 

verging flow, the opposite limits are not approached for a diverging flow; that is, /u^ 

2/2 

does not go to infinity as h>2 zero but approaches a limiting value of 5. This 

occurs in divergii^ flow because transverse as well as longitudinal velocity fluctuations 
are present when the vortex axes are alined transversely, while in convergii^ flow, the 
longitudinal velocity fluctuations approach zero as the vortex axes become alined longi- 
tudinally. Included in the plot (fig. 6) for comparison is the curve for for a sud- 

den contraction without viscosity as obtained from reference 4. 

Relative intensity ratios for turbulent components corrected for viscous decay are 
plotted in figure 7. To obtain the ordinates in this figure, the ordinates for the solid 
curves in figure 3 were corrected to eliminate decay by dividing them by the ordinates for 
the dashed curves, which are for a pure viscous decay. The result, after the square root 
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Ftgure 7. - Relative intensity ratios for turbulent components corrected 
to eliminate decay. 


has been extracted, is divided by U/Uq to give intensity ratios relative to the local mean 
velocity. Curves obtained from reference 4 for a sudden contraction and no viscosity 
are included in the plot in figure 7 and show fair agreement with the present analysis. 


The curve in figure 7 for < 



can be related approxi- 


mately to the heat transfer between a cone wall and a fluid, which occurs, for instance, 
in a cooled rocket nozzle. The comparison is made by using the following argument, 
which is based on the momentum-heat-transfer analogy: Except very close to the wall, 
the turbulent heat transfer is large compared with the molecular heat transfer, so that 
the total radial heat flow per unit area q is approximately pCpC dT/dr. The radial 
temperature gradient dT/dr at a particular radius is assumed proportional to AT/6, 
where 6 is the boundary layer thickness. The eddy diffusivity e is replaced by the 
product of a transverse velocity fluctuation and a mixing length which is assumed propor- 

0I/2 2 ^^^ 

tional to 6, so that e ~ u, 6. Then the heat-transfer coefficient h = q/AT ~ PCpUg 


> 1/2 


/u. 


2 ^ 

or the Stanton number St = h/pUc^ ^ u„ 

P 2 

It is assumed that the change in turbulent intensity along the cone in the boundary 
layer is determined primarily by the normal strain rather than by the shear, as it might 
be when the strain is very rapid. (The initial turbulent intensity in the botmdary layer 
at the entrance of the cone would, of course, be determined by the shear in the upstream 
boundary layer. ) By using the preceding relation for Stanton number, the ordinate curve 
for i = 2 in figure 7 can then be replaced by St/Stg. That curve, which has been 
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a - 30° 1 Experf mental data for flow fn cone, 
a - 60° J Mach number < 0. 2 (ref. 12) 

Present analysis based on normal strain (curve for i = 2 
in fig. 7) 

Calculated convective heat transfer for boundary layer in 

cone with effects of normal strain neglected (eq. (51)) 



Figure 8. - Approximate calculation of heat transfer between cone wall 
and fluid and comparison with experiment. Angle, a degree, in 
figure is angle between cone generator and flow direction. 


corrected to eliminate viscous decay, is used because the inhomogeneous shear which 

occurs in a boundary layer will normally offset the viscous decay. The curve, replotted 

in figure 8, shows that Stanton number decreases rapidly along a converging cone. If, 

by contrast, the Stanton number is calculated from local conditions in a boundary layer 

without considerii^ normal strain, the decrease is much more gradual. In that case the 

0 2 

Stanton number is roughly proportional to (pUD)" (assuming that 5 ~ D), and 

St/St(j = (1 + (51) 

The curve for equation (51) is the dashed curve in figure 8. 

The heat transfer in the boundary layer of a rocket nozzle may, of course, be more 
complicated than the case considered herein, where changes are assumed to be governed 
by the normal strain. The effects of shear and variable properties, as well as of nor- 
mal strain, may not be negligible. However, experimental data for heat transfer in 
cooled conical nozzles (refs. 12 and 13) indicate trends very similar to those obtained 
herein. Data from reference 12 for M < 0. 2 in the conical portions of two nozzles are 
plotted in figure 8. Comparison of the data with the solid and dashed curves seems to 
indicate that changes in Stanton number along the cone are much more dependent on the 
normal strain than on the shear in the boundary layer. When the data were plotted, it 
was assumed that the entrance of the cone corresponds to x* = 0. Although there is 
some uncertainty as to the point in the analysis which corresponds to the entrance of the 
cone in the experiment, it turns out that the results are insensitive to the point chosen. 

H, for instance, x? = -0. 5 were chosen as the starting point instead of x’ = 0, the 
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results would be nearly identical. 


Turbulence Spectra 

2 ^ 2 ^ 

Spectra of and Ug are plotted in figures 9 and 10 and show how contributions 
2 * 2 * 

to Uj and Ug are distributed among dimensionless wave numbers (or reciprocal eddy 
sizes)K*. Plotted in this way the curve for no convergence or divergence (x' = 0) does 
not change with x. Thus, comparison of the curves for various values of x’ with the 
curve for x’ = 0 shows how convergence or divergence affects the spectrum at a given 
position in comparison with the spectrum at the same position with no strain. For con- 

2 * 

vergii^ flow (negative x’) the contributions to u^ occur at smaller wave numbers 

2 * 

(larger eddies) than they would for no convergence, whereas contributions to Ug move 

to higher wave numbers. For diverging flow, by contrast, both and 

to lower wave numbers than they would for no divergence. — 

2 * 

Of some theoretical interest is the extreme asymmetry of the spectrum of u^ for 
negatively large x* . This effect has been observed previously (refs. 6 and 9), but it is 
much more pronounced here, possibly because the strain rate increases sharply as x? 
becomes more negative. As in the previous cases the asymmetry is associated with a 
spectral transfer term that depends on the mean strain rate. That the effect is asso- 
ciated with a transfer term (eq. (48)) was verified by solvii^ the spectral equation with 
the transfer term omitted. The spectrum obtained was found to be nearly symmetrical. 

A plot of the transfer term associated with the mean strain Tjj is given in figure 11. 
The net area under each curve is zero, in agreement with equation (46). The curve for 
negative x’ is predominately negative at low wave numbers and positive at higher ones, 
so that energy is in general transferred from low wave numbers to higher ones. The 
curve also indicates that a small amount of reverse transfer to low wave numbers takes 
place in the very low wave number region. 

The fact that the energy transfer in the longitudinal component of the turbulence is 
primarily from low to high wave numbers is possibly surprising, since according to a 

2 

simplified theory it might be assumed that the Uj^ component is produced by vortices that 
are alined transversely. The axes of those vortices would tend to be shortened in a con- 
verging flow. However, the proportion of vortices alined in the transverse direction is 

p 

small in a converging flow. Most of the contribution to u| probably comes from oblique 
vortices, and in those vortices the energy transfer can be in the direction indicated in 
figure 11. For diverging flow (positive x’) the energy transfer is in the opposite direc- 
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Dimensionless wave number, 

Figure 11. - Transfer spectra x' is an axial distance parameter. 


2.4x10'5 



Dimensionless wave number, «• 

Figure 12. - Energy and dissipation spectra for u^ for x'--0.99. 
(Dissipation spectrum normalized to height of spectrum. ) 


tion, that is, from high to low wave numbers. In that case the effect of the energy trans- 

2 * 

fer on the shape of the Uj spectrum seems to be small. 

For the converging case, figure 11 shows that the positive area in the high wave 
number region is spread out over a wide range of wave numbers. As x* becomes more 
negative this range widens still more. This elongated positive area of energy transfer is 

2 * 

responsible for the long tail on the spectrum of u| . To carry the effect to the extreme, 

2 * 

the energy spectrum of u^ for x’ = -O. 99 is plotted semilogarithmically in figure 12. 
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2 

Included also in the plot is the dissipation spectrum, which is proportional to k rj/.^y 
The energy and dissipation regions in this case show essentially complete separation and 
are separated by a pseudo- inertial subrange in which energy transfer is the dominating 
process. This inertial subrar^e is termed pseudo because it occurs only in one compo- 
nent of the turbulence and because it is produced by inertial effects associated with the 
mean strain rather than with the triple correlations that are usually considered to be 
responsible for an inertial subrange (ref. 8). Figure 12 is, however, a good illustration 
of a case in which essentially complete separation of energy and dissipation regions 

occurs. — ^ 

2 * 

Figure 13 shows a log-log plot of the spectrum of u^ for x* = -0. 99. In this plot 
is proportional to k*~^ over about four decades of k*. The region in which the 
curve begins to fall off more rapidly than k*~^ is roughly the region in which dissipation 
effects become important. For x’ = -1, the dissipation region would be moved to infinity 
and x}/* would vary as k*~^ over the entire range of wave numbers. The present results 
for a weak turbulence with large mean strain rates differ from the usual Kolmogorof -5/3 
power spectrum that appears to apply at very high Reynolds numbers (ref. 14). As shown 
by the present results, however, a -5/3 power region in an energy spectrum is not nec- 
essary for complete separation of the energy and dissipation regions. In fact, any power 
between 0 and -2 will do as well, since the dissipation spectrum is obtained by multi- 
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2 

plying the energy spectrum by k . 

Vorticity spectra were also calculated, and representative results for spectra of 

2 * 

(x >2 are plotted in figure 14. As x* becomes more negative, the spectra move to higher 

2 * 

dimensionless wave numbers. Contributions to cOg also become spread out over a much 
wider rai^e of wave numbers or vortex sizes. As x’ increases positively, the spectra 

2 * 

move to lower dimensionless wave numbers. The trends for the spectra of Wj (not 
shown) are similar to those shown in f%ure 14 with the exception that the shapes of the 
spectra do not change appreciably with strain. Thus, the vortices in converging flow 
tend to be smaller at a given x than they would be for no convergence, whereas they 
tend to be comparatively large for diverging flow. 


SUMMARY OF RESULTS 

Near the virtual origin of the turbulence, for both converging and diverging flow 
through a cone, all turbulence components decreased along the flow because of viscous 
dissipation. For a converging flow the effect of the distortion was to tend to aline the 
turbulent vorticity in the direction of mean flow. This caused the longitudinal component 
of the velocity variance to decrease toward zero and the lateral components to increase. 
For a diverging flow the vorticity tended to be alined in the transverse direction. In that 
case the ratio of longitudinal to transverse velocity variance ultimately approached 5. 
When the results for turbulent intensity were corrected to eliminate viscous decay and 
divided by local mean velocity, both longitudinal and transverse components decreased 
along a converging flow and increased along a diverging flow. The results were related 
approximately to heat transfer between a cone wall and a fluid and gave trends very 
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similar to those observed experimentally. 

Turbulent vorticity spectra showed that the turbulent vortices in a converging flow 
tend to be smaller at a given position than they would be for no convergence, whereas 
those in a diverging flow tend to be relatively large. A mean-strain transfer term in the 
spectral equation for the longitudinal component of the turbulence transferred energy in 
the high wave number direction for a converging flow. This transfer caused the spectrum 
for the longitudinal component of the energy to become strongly asymmetric. Near the 
cone apex essentially complete separation of the longitudinal energy and dissipation 
spectra was obtained. A power spectrum of -1 rather than -5/3 that extended over a 
considerable range of wave numbers was also noted for the longitudinal energy component. 
Thus, a -5/3 power spectrum is not a necessary consequence of the separation of the 
energy and dissipation regions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 9, 1966, 

129-01-09-07-22. 
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APPENDIX - SYMBOLS 


B 

defined by eq. (41) 

U.U.* 
1 ] 

defined by eq. (37) 

D 

local diameter 

X 

defined by eq. (21) 

g 

defined by eq. (39) 

X* 

defined by eq. (28) 

H 

defined by eq. (40) 


(x - Xq) tan 0 ! 

h 

heat transfer coefficient 

A. 


n 

defined by eq. (38) 

Xi,(Xi)m 

defined in fig. 1 

*^0 

constant that depends on 

x,(xi)^ 

longitudinal coordinate 

initial conditions 

a 

angle between cone gener 

M 

local mach number 


ator and flow direction 

p 

instantaneous pressure 


(fig. 2) 

r 

radius 


defined by eq. (45) 


defined in fig. 1 

6-. 

Krone cker delta 

r 

1] 

Tn 

distance from centerline 

e,cp 

angular coordinates (see 

0 

to cone generator at Xq 


eq. (31)) 


(fig- 2) 


wave number component 

S,S(1) 

aUi _ _ 3(Dj)^ 

K* 

defined by eq. (29) 

axj ax 


defined by eq. (43) 

®(2)’ ®(3) 



- Xq)^ 

^ij u^j 




defined by eq. (8) 



defined by eq. (9) 

Til 

longitudinal transfer term 

V 

kinematic viscosity 

defined by eq. (48) 

P 

density 

t 

time 

<Pij 

defined by eq. (7) 


mean velocity component 

(P*j 

defined by eq. (30) 

Ui,(Ui)^,U 

mean longitudinal velocity 


defined by eq. (33) 


component 



fluctuating velocity 
component 
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IJ 






CO-0). 


defined by eq. (42) 
turbulent vorticity variance 
.V2(,.^)7/2 


CO-CO- 
1 ] 


'’o 


Subscripts: 

a conditions upstream of contrac- 

tion at point where mean veloc- 
ity begins to vaiy (fig. 4) 


w without effects of strain 

0 at virtual origin of turbulence 

where turbulent energy would 
be infinite; it is assumed that 
turbulence is isotropic at Xq 
and that strain begins to act 
there 

1 in longitudinal direction 

2,3 in transverse direction 

Superscripts: 

' at point P (see fig. 1), except 

for dimensionless x* 


c corrected to eliminate viscous 

decay 

m at point (see fig. 1) 


dimensionless quantity 
averaged quantity 
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